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Abstract A new version of improved quantum molecular dynamics model that includes standard Skyrme 
interactions has been developed. Based on the new code, four commonly used parameter sets, SLy4, SkI2, SkM“ and 
Gs are adopted in the improved quantum molecular dynamics model and the isospin sensitive observables, namely 
isospin transport ratios, single and double ratios of the yields of neutrons and protons are investigated. The isospin 
transport ratios are strongly sensitive to the slope of symmetry energy, and are not very sensitive to the nucleon 
effective mass splitting. On the other hand, the high energy neutrons and protons yields ratios from reactions at 
different incident energies provide a good observable to the momentum dependence of nucleon effective mass 
splitting. By comparing our calculations with the data, we find that the constrained L value (the slope of density 
dependence of symmetry energy) is about ~46 MeV when the Skyrme type interaction is considered in transport 


models, and the isospin diffusion data prefer to m > m, , but it is not a strong constraint with deep y? minimum. 
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1 Introduction 


Knowledge about the symmetry energy which is 
defined as the difference in the binding energy 
between the pure neutron matter and symmetric 
nuclear matter is important for understanding not only 
the nuclear structure and nuclear reactions but also 
many critical issues in astrophysics. Recent 
observables used to constrain the symmetry energy 
range from isospin diffusions in heavy ion collisions 
(HIO, to experiments that measure nuclear 
properties such as neutron skin, Pygmy Dipole 
Resonance, masses of Isobaric Analog States, and 
nuclei masses in Finite Range Droplet Model!!>4), 
Until now, a general consensus on the symmetry 
energy, So, at saturation density and its slope, 
L (L=3p,d8(p)/dp VER" where S(p) is the density 
dependence term of the symmetry energy), were 
obtained. But, the results still have large uncertainty 


and it needs further constraints. 


The constraints from heavy ion collisions are 
obtained by comparing the data of isospin sensitive 
observables with the predictions from transport models 
which include the density dependence of symmetry 
energy as an input variable. Tighter constraints can be 
obtained from improved experiments or from 
improved theoretical models by including the missing 
physics. 

In addition to the dependence on nuclear 
density, the symmetry potential also depends on the 
momentum of the nucleons in the system. Such 
dependence leads to the splitting of the nucleon 
effective mass, especially in the high density regions. 
An interesting new question arises as to whether the 
effective mass for neutrons (m, ) is higher than, equal 
to, or lower than that for protons (m, ), and how does 
the effective mass splitting changes with the 
momentum of nucleon. In fact, our lack of reliable 
knowledge about the neutron-proton effective mass 
splitting has been hindering the accurate extraction of 


S(p) from nuclear reactions by using transport models, 
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because the compressed nuclear system is excited due 


to violent nucleon-nucleon collisions and the 
contributions from the momentum dependence of 
symmetry potential cannot be entirely described with 


[1.25,26] 


local density approximation There are some 


efforts on constraining the nucleon effective mass 
splitting by analyzing nuclei optical potential datal!?7) 
and the m,<m, has been ruled out around normal 
density and Fermi momentum. Furthermore, it would 
be interesting and important to give constraints on the 
momentum dependence of effective mass splitting 
from heavy ion collisions data. In this paper, we 
investigate such dependence in a transport model by 
choosing different Skyrme effective nucleon-nucleon 
interactions. 


2 New version of ImQMD 


The Quantum Molecular Dynamics Model (QMD) 


represents the individual nucleons as Gaussian 


"wave-packet" with mean values that move in 
according the Ehrenfest theorem; i.e. Hamilton's 


equations"! 


. At China Institute of Atomic Energy 
(CIAE), we developed and successfully applied a new 
QMD code, labeled IMQMD which includes mean 
field potentials calculated using a Skyrme energy 
density function with options for different forms of the 
density dependence of the symmetry potential, to 
study heavy ion reactions ranging from Coulomb 
barrier to 400 MeV per nucleon'*!, With these 
modifications, ImQMD has successfully described the 
multiplicity of reaction products, collective flows and 
stopping powers in intermediate energy Heavy Ion 
Collisions (HICs), but this version lacked to consider 
the nucleon effective mass splitting. 

We improve the mean field part of the 
Improved Quantum Molecular Dynamics (ImQMD05) 
code by including the isospin dependent Skyrme like 
momentum dependence interaction as the form 
Ung © [apap f (rp) f(np'\(p-p’) - The form of 
Uma iS 


Una > Und (PT) +U yg (PTa) + Una (2,7, ) 
=C,[d’ pd’ p'f (r.p) f(r. P')(p-P') 
+D, |d pd’ p'lf, (r.p) f, (rP) PP’) 


+f, (rp) f,(r,P')(p-P’) 1 


(1) 


where f is the nucleon phase space distribution 


function: 


er aof Pool en) (3) 
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The coefficients Cy and Do can be determined with 


following relationship: 


1 

Cy = Tegal (24%) +h (2+0)1 (5) 
1 

D, =zzzlh (2% +1)-5 (24 +1) © 


In the new version of IMQMD code (ImQMD-Sky), 
the potential energy U is 


U = Uoc + U ma T U sou z (7) 
The nuclear contributions are represented in a form 
with Uema = Wiocmad r. The energy density toc is 


2 ntl 
_&P 4 B P- m Esur [Ve] 
22 Ptl py 2P 
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Here, 6 is the isospin asymmetry. &=( Pa, -pp)/ (Put Pp ), 
fn and p, are the neutron and proton densities, 
respectively. The coefficients of œ, 2, sur 8surisos Asym 
and Bym can be obtained by the standard Skyrme 
interaction parameters as in previous work), Ueout 18 
Coulomb energy. The mean fields acting on these 
wavepackets are derived from potential energy. 

These calculations use isospin-dependent 
in-medium nucleon-nucleon scattering cross sections 
in the collision term and Pauli blocking effects as 
described in Ref.[28]. With this new version of 
Improved Quantum Molecular Dynamics code, we can 
directly test the validity of Skyrme interaction 
parameters which have been wildly used in nuclear 
structure studies in heavy ion collisions, and provide 
constraints on the density dependence of symmetry 
energy and n/p effective mass splitting. 


Based on the Skyrme interaction, one can get 
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the density dependence of symmetry energy, neutron 
and proton effective mass and the symmetry potential 
for cold nuclear matter on the mean-field level as 


follows: 
TE. 7 2 
Sup ep 
(mi /m\ =[1+4(C,p+Dyp, )m], (10) 
U,-U 
U m = e > =2A mP +2B mP” +2DmpE,.(11) 


Eq.(9) is the density dependence of symmetry energy 
for cold nuclear matter, the first term comes from the 
contribution of kinetic energy part, the second and 
third terms are from the two-body and three body 
terms in Skyrme interactions. The last term in Eq.(9) 
comes from the momentum dependent interaction term 
in Skyrme interaction, and C,,m is determined as in 
Ref.[28]. Eq.(10) is the effective mass where p;is the 
neutron or proton density for t=n, p. Eq.(11) is the 
symmetry potential also known as the Lane potential 
which gives the strength of the symmetry potential in 


asymmetric nuclear matter. 
3 Results and discussion 


In the following studies, we choose four Skyrme 
interaction parameter sets, SLy4, SkI2, SkM” and Gs 
which have similar effective mass and 
incompressibility’”?*"!, i.e., m'~0.740.1, Ko=230+20 
MeV, and have been widely used in the studies of 
nuclear structure and neutron star. The SLy4 and SkI2 
were fitted to properties of neutron matter, neutron star 
and the ground-state variables of neutron-rich heavy 
nuclei, and the neutron effective mass is less than the 
proton effective mass, m; <m. We choose these two 
interactions because the slopes of symmetry energy L 
are very different; L is 46 MeV for SLy4, and 
104 MeV for SkI2. The SkM’ and Gs were fitted to 
binding energies of finite nuclei and actinide fission 
barriers, and the neutron effective mass is larger than 
the proton effective mass, m >m. L for SKM’ is 
46 MeV, and 93 MeV for Gs. By analyzing the results 
from calculations using SLy4, SkI2, SkM” and Gs and 
comparing to data, we hope to understand the 


sensitivities of the interactions on the density 
dependence of symmetry energy and neutron proton 
effective mass splitting simultaneously and get their 
constraints. In this paper, we focus our discussion on 
the '!*!4Sn4!!?!748n reactions at 50A MeV which 
have been measured in the MSU/NSCL group. 

Figure 1 shows the density dependence of 
symmetry energy (left panel) and the energy 
dependence of the Lane potential (right panel), for 
cold nuclear matter. At subsaturation density, the 
strengths of symmetry energy obtained with SLy4 (top 
solid lines) and SkM“ (bottom solid line) are stronger 
than that obtained with SkI2 (top dashed lines) and Gs 
(bottom dashed lines). Separation of the two groups of 
(solid and dashed) lines is strongly correlated with the 
L values. Smaller L values yield higher symmetry 
energy at subsaturaion densities while the opposite is 
true at the suprasaturation density regions. The right 
panels of Fig.1 show the Lane potentials for the above 
four interactions at 0.5, (top panel) and at p (bottom 
panel) as a function of nucleon kinetic energy. The 
energy dependence of symmetry potential for SLy4 
and SkI2 are positive and increase with the energy of 
nucleons. However, those for SkM” and Gs decrease 
with the kinetic energy and become negative above 
200 MeV. When m, <m 
SkI2, the neutrons feel stronger repulsion than protons 


p» as in the case of SLy4 and 
and the effects increase with the kinetic energy 
increasing. The opposite is true for SkM” and Gs 
where m, >m, the protons feel stronger repulsion 
with the kinetic energy increasing. 

The symmetry potential (Lane potential) also 
gives an accurate estimate for the difference of the 
force between neutron and proton feeling in 
asymmetric nuclear matter, and directly influences the 
n/p ratios of emitted nucleons. The larger the Lane 
potential is, the larger the n/p ratio is. At high energy, 
one can expect that the n/p ratios with m, <m, are 
greater than that for m; >m, because the Lane 
potential for m, <m, is larger than that for m >m, 
cases. Furthermore, we can see that a crossing of two 
symmetry potentials for SLy4 and SkM“ appears at 
Es ~ 30 MeV for p=0.5p9 and Eps ~ 70 MeV 
for p=. Below the crossing energy, the Lane 
potential for m; <m; cases is less than that for 


m >m, , we can expect that the n/p ratios with 
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m,<m, are less than that for m,>m, if their 


differences on the Lane potential are large enough. 


P! Py E I MeV 


Fig.1 (Color online) (a) Density dependence of symmetry 
energy, (b) Energy dependence of the Lane potential at p=0.5 
and (c) p=, for the parameters SLy4 (red solid lines), SkI2 
(red dashed lines), SkM (blue solid lines), Gs (blue dashed 
lines). 


Before analyzing the isospin sensitive 
observables, we firstly check the calculated results on 
charge distribution for 14Sn+!Sn with ImQMD-Sky. 
4en4'4Sn at 
Eveam=0A MeV obtained with four Skyrme parameter 
sets, SLy4, SkI2, SkM and Gs. It is clear that the 


charge distribution, which is considered as an isoscalar 


Fig.2 shows the charge distributions for 


observable, is weakly sensitive to the slope of 
symmetry energy and nucleon effective mass splitting 
among the selected parameters. 


e 
124Sn+124Sn, 50 AMeV, b=6 fm 


Yield / Arb. 


he 
*Mennoeeeeeseagy, 


Fig.2 (Color online) Charge distributions for '“Sn+'™Sn at 
E/A=50 MeV and b=2 fm for four Skyrme parameter sets, SLy4, 
SkI2, SkM” and Gs. 


studies 68211 the 


In previous density 
dependence of symmetry energy has been probed 
using peripheral collisions to measure the isospin 
diffusion and isospin transport ratios as a function of 
rapidity. Moreover, even though there are data from 
n/p yield ratios and flow, the constraints on the 
symmetry energy from heavy ion collisions rely 


heavily on isospin diffusion data!®!", 


In this paper, we simulated the collisions of 
Menon P4Sn+! Sn, !?Sn+Sn, and ony! gn 
reactions at beam energy of 50 MeV per nucleon using 
the ImQMD-Sky code. 64,000 events are obtained for 
each reaction at b=6 fm. In the left panel of Fig.3, we 
plot the isospin diffusion transport ratios obtained with 
SLy4, SkI2, SkM” and Gs interactions. As in previous 


studies!®®! 


, we analyze the amount of isospin diffusion 
by constructing a tracer from the isospin asymmetry of 
emitting source which includes all emitted nucleons (N) 
and fragments (frag) with the 


cut (v~s >0.5ve™ ). The shaded regions are the 


z beam 


velocity 


isospin diffusion data obtained by constructing the 
isospin transport ratios with isoscaling parameter X=a, 
near the projectile rapidity regions. Faster 
equilibration occurs for smaller L values which 
correspond to larger symmetry energies at 
subsaturation density. Our results show that the R; 
values for SLy4 (solid circle) and SkM“ (solid squares), 
both with L=46 MeV, are smaller than those for SkI2 
(open circle, L=104 MeV) and Gs (open square, L=93 
MeV) with similar effective mass splitting. The overall 
effect of mass splitting on isospin diffusion is small. 
When m,<m,, the isospin diffusion process is 
accelerated at subsaturation densities due to the 
stronger Lane potential. R; values increase in the order 
of R(Sly4: L=46 MeV, m, <m )< R(SkM’: L=46 
MeV, m > m, )< R(SkI2: L=104 MeV, m < m, )< 
R{Gs: L=93 MeV, m,>m,). By comparing to the 
data, we found that the R; obtained with both SLy4 and 
SkM” (L=46 MeV) fall within the experimental 
uncertainties among the interactions we adopted. The 
data prefer the softer symmetry potentials with lower L 
values. 

We also compare the results of the calculations 
to R; as a function of the scaled rapidity y / ygn, with 
SLy4, SkI2, SkM” and Gs interactions as shown in the 
right panel of Fig.3. The star symbols in the right 
panel are experimental data obtained in Ref.[13] for 
peripheral collisions. This transport ratio was 
generated using the isospin tracer X=In[Y(Li)/ Y(‘Be)], 
where YCLi)/Y( Be) is the yield ratio of the mirror 
nuclei, Li and Bel’, For comparison, the 
ImQMD-Sky calculations of R(X=6y,frag) are plotted 
as lines for b=6 fm. The interactions with smaller L 


values, SLy4 and SkM” (solid lines) agree with the 
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data better especially in the high rapidity region. 
However, y’ analysis suggests that the data fit the 
results from SkM” interactions with nč > m’, 
effective mass splitting. Current study show that the 
sensitivity of the isospin diffusion to slope parameter 
is consistent with previous studies but the sensitivity to 
effective mass splitting is not large. 

0.8 


0.6 


0.2 


0 Sly4 Ski2 SkM* Gs 0 0.5 1.0 
Yl ee 


Fig.3 (Color online) Left panel: The isospin diffusion 
transport ratios obtained with SLy4, SkI2, SkM” and Gs. The 
shadow region is the data from Ref.[3]. Right panel: The 
calculated results on isospin transport ratios as a function of 
rapidity for SLy4, SkI2, SkM“ and Gs. The star symbols are the 
data from Ref.[13]. 


Rizzo®! has shown that the neutron/proton 
yield ratio, Y(n)/Y(p), as a function of p, from central 
collisions is a robust observable to study nucleon 
effective mass splitting. In Fig.4, we plot the 
R(nlp)=Y(n)/Y(p) ratio for 12Sn+!!°Sn (left panel) and 
124, | 124 

Sn+ 
gate 70°<6.m<110°. The lines connecting the solid 


and open circles correspond to m, < m, case, and the 


Sn (middle panel) at b=2 fm with angular 


lines connecting the squares correspond to m, >m, 
case. Not surprisingly, the Y(n)/Y(p) ratios are larger 
for the neutron rich system, 14Sn+!”*Sn, in the middle 
panel. Consistent with Refs.[26,27], the differences in 
the Y(n)/Y(p) ratios between the m,<m, (circles) 
and m, > m, (squares) increase with kinetic energy of 
the nucleons. In this high energy region, the single 
ratios (left and middle panels) are more sensitive to the 
effective mass splitting than the isospin diffusion 
observables. This is because the contributions from 
momentum dependent part of symmetry potential play 
more and more important roles on the Y(n)/Y(p) ratios 
when the relative momentum of nucleons increase. As 
shown in the right panel of Fig.1, the Lane potentials 
become larger for m,<m, at high kinetic energy, 
and cause enhanced neutron emissions. Thus, the 


information about the momentum dependence of 
nucleon effective mass splitting is contained in the 
energy spectral of pre-equilibrium transverse emission 


of nucleons at earlier stage of reactions. 


tr] 2.5 
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Fig.4 (Color online) Left panel is the calculated results for 
Y(n)/Y(p) as a function of kinetic energy for 112Sn+!!?°Sn at b=2 
fm with angular cuts 70°<@,,,<110°; middle panel is the 
Y(n)/Y(p) for 246 n4)4Sn. Right panel is the calculated results 
for DR(n/p) ratios as a function of kinetic energy. The 
calculated results are for SLy4 (red solid lines), SkI2 (red 
dashed lines), SkM” (blue solid lines), Gs (blue dashed lines) 
and NRARP (green lines). Star symbols are the data", 


The calculations with SLy4 interactions agree 
with the double ratios data [4], DR(n/p)=(Ya(n)V/Y4(——)) 
(Yp,(n\/Yp(p)), where A='“Sn+'“Sn and B=''?Sn+ 
12S, from Ref.[4] 
uncertainties especially for ratios of the high energy 


which have very large 


nucleons. Thus, exact constraints on effective mass 


splitting from heavy ion collisions need 
re-measurements of the data with high quality. Since 
the effect of mass splitting should be larger at higher 
densities, most likely, one needs data on the DR(n/p) 
ratios at different incident energies which should give 
constraints on the momentum 


more stringent 


dependence of n/p effective mass splitting by 


comparing with improved transport models 


calculations. 


4 Conclusion 


In summary, we have used a new version of the 
improved quantum molecular dynamics code, which 
can accommodate Skyrme interaction parameters, to 
study the isospin sensitive observables, such as isospin 
diffusion, isospin transport ratios as a function of 
rapidity, single and double neutron proton yield ratios. 


We find that the neutron proton effective mass 
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splitting plays an important role on the n/p ratios of 
transverse emitted nucleons at high kinetic energy. 
The mass splitting affects isospin diffusion, but the 
effects are not strong. Our results show that the 
constrained L~46 MeV value. The isospin diffusion 
data prefer to the neutron effective mass is greater than 
proton effective mass, but it is not a strong constraint 
with deep x minimum. We also show that the single 
and double n/p ratios from high nucleon energies are 
more sensitive to the mass splitting effects. The results 
also suggest that the effective nucleon masses should 
be further studied at different incident energies. 
Finally new neutron and proton spectral data with 
much smaller uncertainties than previous data at 
different beam energies may allow one to determine 
the sign and the momentum dependence of the mass 
splitting. 
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